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EFFECT OF SPANWISE VARIATION OF TURBULENCE ON 


THE NORMAL ACCELERATION OF AIRPLANES WITH 
SMAIL SPAN RELATIVE TO TURBULENCE SCALE 


By Kermit C. Pratt 


ABSTRACT 


A rigid airplane with an unswept wing is analyzed. The results show 
that the power spectrum, relative to that for a one-dimensional turbulence 
field, is significantly attenuated at the higher frequencies even for air- 
planes with arbitrarily small ratios of span to scale of turbulence. This 
attenuation is described by a simple weighting function of frequency that 
depends only on aspect ratio. The weighting function, together with the 
attenuation due to the unsteady flow of gust penetration, allows the determi- 
nation of the average rate of zero crossings for airplanes having very small 
spans without recourse to an integral truncation that is often required in 
calculations based on a one-dimensional turbulence field. 


SUMMARY 


The normal acceleration response of a rigid airplane, with an unswept 
wing, to a two-dimensional turbulence field is analyzed. The analysis 
features an atmospheric turbulence spectrum in terms of wave numbers in the 
direction of flight and along the span. The spanwise distribution of aero- 
dynamic forces is approximated by strip theory. Spectral integrations with 
respect to spanwise wave numbers are obtained in closed form for asymptotic 
conditions attending small and large xrave numbers in the direction of flight. 
This approach reveals some significant and fundamental effects of the two- 
dimensional turbulence field and parametric relationships for airplanes with 
small span-to-turbulence-scale ratios (in the range from about thirty percent 
to vanishingly small). Results relative to those for a one-dimensional 
turbulence field show that the power spectrum is attenuated at the higher 
frequencies even for airplanes with arbitrarily small ratios of span-to- 
turbulence scale. The effect is accounted for by a simple weighting function 
of frequency that depends only on aspect ratio. This attenuation, together 
with that from the unsteady flow of gust penetration, permits the determination 
of the average number of zero crossings without recourse to an integral 
truncation that is often required in one -dimensional turbulence response 
calculations. 
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EFFECT OF SIANVISE VARIATION Of TUR3ULENCE ON 
THE NORMAL ACCELERATION OF AIRPLANES WITH 
SMALL SPAN RELATIVE TO TURBULENCE SCALE 

By Kermit G. Pratt 


ABSTRACT 


A rigid airplane with an unsvept wing is analyzed. The results show 
that the power spectnun, relative to that for a one- dimensional turbulence 
field, is significantly attenuated at the higher frequencies even for air- 
planes with arbitrarily 3mall ratios of spar, to scale of turbulence. This 
attenuation is described by a simple weighting function of frequency that 
depends only on aspect ratio. The weighting function, together with the 
attenuation due to the unsteady flow of gust penetration, allows the determi- 
nation of the average rate of zero crossings for airplanes having very small 
spans without, recourse to an integral truncation that is often required in 
calculations based on a one-dimensional turbulence field. 


SUMMARY 


The normal acceleration response of a rigid airplane, with an unswept 
wing, to a two-dimensional turbulence field is analyzed. The analysis 
features an atmospheric turbulence spectrum in terms of wave numbers in the 
direction of flight and along the span. The spar.wise distribution of aero- 
dynamic forces is approximated by strip theory. Spectral integrations with 
respect to spanwise wave numbers are obtained in closed form for asymptotic 
conditions attending small and large "ave numbers in the direction of flight. 
This approach reveals some significant and fundamental effects of the two- 
dimensional turbulence field and parametric relationships for airplanes with 
small span-to-turbulence-scale ratios (in the range from about thirty percent 
to vanishingly small). Results relative to those for a one-dimensional 
turbulence field show that the power spectrum is attenuated at the higher 
frequencies even for airplanes with arbitrarily small ratios of span-to- 
turbulence scale. The effect is accounted for by a simple weighting function 
of frequency that depends only on aspect ratio. This attenuation, together 
with that from the unsteady flow of gust penetration, permits the determination 
of the average number of zero crossings without recourse to an integral 
truncation that is often required in one-dimensional turbulence response 
calculations . 


INTRODUCTION 


The methodology for determining the theoretical statistical averages 
that describe the symmetric responses of an airplane flying in continuous 
turbulence i»as originally developed on the basis of a one-dimensional 
turbulence model*, i.e., turbulence velocities that vary in the direction 
of flight and are invariant a?.ong the span (See bibliography in reference l). 
The procedure is currently in wide use and provides useful information for 
a large number of problems. There is, however, an inherent awkwardness in 
that in some instances (e.g., for an essentially rigid airplane) the calcula- 
tion of the average rate of zero crossing requires the truncation of the 
upper limit of an integral (ref. 2). This injects the effect of engineering 
judgment and so causes some uncertainty in the results. 

Methods for determining the effects of spanwise variations in turbulence 
(two-dimensional turbulence) have been developing for many years. The 
earliest investigators, references 3-9, were primarily concerned with 
the mathematical concepts for treating the two-dimensional problem. Attention 
was given mostly to mean-square responses; particularly of the lift due 
directly to the turbulence and of normal acceleration to a lesser degree. 

In view of the then newly-acquired ability, it is not surprising that 
emphasis was placed on fairly large values of the ratio of wing span to the 
turbulence scale (up to unity and beyond). There was a general tendency to 
expect the results from the two-dimensional turbulence model to approach 
those for the one-dimensional model as the ratio of span to turbulence scale 
approaches zero. 

In the past few years, however, there has been a growing appreciation of 
the significant effects of the two-dimensional turbulence field on the 
symmetric responses of airplanes having spans that are less than the order 
of ten percent of the turbulence scale. (This figure is representative of 
the great majority of airplanes at cruising altitudes if a value of 762 meters 
(2500 feet) is accepted.) In reference 10, spectra of normal acceleration 
and symmetrical wing bending moment for an airplane having a span-to-scale 
ratio of five percent were shown to be attenuated at the higher frequencies 
by the two-dimensional turbulence model relative to those for the one- 
dimensional model. The attenuation was to such an extent that when combined 
with the unsteady effects of gust penetration, the values of the average rate 
of zero cr ing could be calculated without recourse to the truncation of 
integrals required by the one-dimensional model. The mean square values of 
the responses were only slightly smaller than those from the one-dimensional 
model. At about the same time similar results and conclusions were reached 
independently in the work reported in reference 11. Also, the effects of 
the two-dimensional turbulence model on the acceleration of the rigid body 
• and elastic modes of the Concorde supersonic transport were examined in 
reference 12. A more recent paper, reference 13, emphasizes the importance 
of the role of two-dimensional turbulence in calculations of the average 
rate of zero crossings. Presumably, there are many others in the world-wide 
aerospace industry whose unpublished data from computer calculations have 
made them aware of these effects. Nevertheless, there are still numerous 
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engineers who are not familiar with the differences in effects of one- and 
two-dimensional turbulence models. One of the pui’poses of this paper is to 
further the dissemination of knowledge of these differences. 


The principal purpose* however, is to describe some insights into the way- 
in which the two-dimensional turbulence model generates the differences from 
the one-dimensional model and the method by which these insights were 
obtained. Basic relationships are frequently lost in an avalanche of numbers 
generated by a complicated computer program. Such results are usually com- 
puted with high accuracy and are desirable for good design (e.g., ref. l4), 
but generalized effects of changes in airplane or turbulence parameters are 
often obscured. 

A rigid airplane with an unswept wing having a constant chord is analyzed 
The results obtained show that the aspect ratio of the wing, which often 
has not been explicitly expressed in past work, is the single source of the 
significant difference between the results of the one- and two-dimensional 
models for the relatively small span-to-turbulence-scale ratios. It is 
shown too that these differences exist for arbitrarily small values of the 
span-to-seale ratio. The differences are accounted for by a single weighting 
function of frequency, in contrast to the families of functions presented in 
references 11 and 13. Finally, the results indicate that two-dimensional 
gust effects may significantly affect the estimated average rate of zero 
crossings even for very small values of span-to-scale ratio. 
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dimensionless wave number in direction of span, 

Integral scale of turbulence 

normal acceleration increment, g units 

average frequency of zero crossing with positive slope 

weighting function for two-dimensional turbulence 
wing area 
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(k • -S .) 2 

U 2 2L 
time 

airspeed 
gust velocity 


w 


modulus of sinusoidal gust velocity 
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y distance along span 

z normal acceleration 


y* y/(V2) 


K 


g 


p 


a 
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angle of attack due to gust 

mass parameter 

pgScC- e 

L a 

mass ratio 

density of the atmosphere 
root mean square 

kg') two-dimensional power spectrum of turbulence, velocity- 
squared 


\<V> 
K <Y> 

®e 


power spectrum of turbulence, velocity-squared 
effective power spectrum of turbulence, velocity-squared 


^An^l*) power spectrum of normal acceleration increment 


(jJ^ circular frequency in direction of flight 

wave number in direction of flight, radians/unit distance 
fig wave number in direction of span, radians/unit distance 

Subscripts 

An incremental normal acceleration 

w gust velocity 
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breakpoint of asymptotic lines 



APPROACH 


The present approach was originated to check the results obtained during 
the work reported in reference 10. These results, which were obtained from 
a digital computer, indicated that, in general, the mean-square values of 
the responses of the example airplane in two-dimensional turbulence did not 
approach those for the one-dimensional case as the span-to-sc&le ratio 

2 ^ approached zero. An analysis of the phenomenon by the present method 

indicated that, indeed, this characteristic was a reasonable result. 

The approach utilizes an atmospheric turbulence spectrum in terms of 
wave numbers in the direction of flight and along the span, as featured in 
reference 3. This spectrum is more revealing of the response characteristics 
than the averaged and transformed cross-correlation functions or the cross- 
spectra used in references 4-7. The present procedure also utilizes 
simple strip theory for approximating the spanwise distribution of lift due 
directly to the vertical turbulence velocities. There is evidence in 
reference 15 that strip theory provides spanwise distributions of gust 
forces with acceptable accuracy. The combination of this form of turbulence 
spectrum and strip theory, together with the constraint that the ratio of 
wing span to turbulence scale be small (b/L « 1.0) permits spectral 
integrations with respect to spanwise wave numbers to be obtained in closed 
form for asymptotic conditions attending small and large wave numbers in the 
direction of flight. The closed form results provide identification of the 
turbulence and airplane parameters that are involved and their relationships. 

This procedure is applied to a study of the differences noted in the 
past between statistical averages for normal accelerations calculated on 
the basis of one- and two-dimensional turbulence fields for airplanes with 
wing spans that are a small percentage of the turbulence scale. The 
statistical averages involved are the gust normal acceleration sensitivity. 



and the average frequency of zero crossing with positive slope 



The spectrum of normal acceleration, $^(k^') is common to both quantities 
and is the subject of the following analysis. 


ANALYSIS 


of the 


General Relationships The acceleration spectrum is expressed in terms 
e dual -wave-number atmospheric turbulence spectrum, ^wg^l ' ,kg' ) , as 



w g 

where (k^ 1 , kg') .is the acceleration frequency response function. This 

function can be expressed as the product 
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where the function H. is the acceleration due to a sinusoidal lift 

An v 

v e 


coefficient C L and is a function of k^' only. The function H c 

Lg 

is the sinusoidal lift coefficient due solely to the two-dimensional field 
of sinusoidal gust velocities 


^ rr ^ e- 
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( 5 ) 


The spanwise distribution of the aerodynamic lift coefficient due 
directly to the gust field is approximated by strip theory. The 
distribution of gust lift coefficient for an un swept wing having a constant 
chord is 



where 


C ^(o) ^steady lift function due to gust penetration and wake 

Lg 

effects and where from equation (5) 

— // '£t- 

M/J- -r 

V 


L (t cmJ'' z l a * s/*j? I y.) 

ir S • 
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The imaginary term in the parentheses in equation (7) is antisymmetric 
with respect to the wing centerline and so does not contribute to the total 
lift . 


The substitution of equation (7) in (6) yields 

^ fV- r fl* 

* Z Cffa) v 
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The total gust lift coefficient is then 
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With 




(10) 
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the modulus squared for equation (10) is 
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( 11 ) 


The expression for the spectrum of airplane acceleration is obtained by 
substituting equation (ll) in equation (4), and substituting in turn equation 
(4) in equation (3). 



The integral 










is an effective turbulence spectrum. 

The results from the two-dimensional gust field are related to those for 
the one-dimensional gust field by 



where R(k. 1 , 5 =-) is a weighting function which contains the effects of the 

JL dh <y 

4> 

two-dimensional turbulence, and (k^ • ) is the turbulence spectrum 

a wg 

for the one-dimensional field. In terras of the spectrum for the two- 
dimennion turbulence: 
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The combination of equations (13), (l4), and (15) yields the following 
expression for R , 


S? = 


W 7 7JT£V * 


where the integrands in equation (l 6 ) are identified by 


where*, of course 






(19) 










In the following analysis, the characteristics of the weighting 
function R are determined from an examination of these integrand^ 1^ 

and Igt over the range of values of k^' . The Dryden turbulence spectrum 

is used for illustrative purposes. The von Karman Bpectrum can be utilized 
equally well. Results from both spectra are presented subsequently. The 
Dryden spectrum for two-dimensional turbulence (from ref. 10) is 



Lot; i VMuencies .- Consider first the frequency range 
ijwo components of integrand Ig(viz., I the spectra for 


k^' < 1.0. 

= 0 , . 5 , 


The 

and 


IvO, and the function 3") are shown individually as log-log plots in figure 1 
for a value of b/2L arbitrarily chosen to be n/20. The low and high 
frequency asymptotic lines are indicated for all functions. The salient 
feature of the functions is the break frequency between the asymptotes. 


The break frequency for S is 
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the value of S is nearly unity. The maximum Vreak frequency for the 
spectra 3^ is 

-4' ar = z. f/t ■■ /■v? A! -- '•*) 

If k « « kp r _ 

B,I B,S 

or 1/(8 « “ 

D 

or b/2L « .5 6 
then S 1 . 0 
and 

Thus the one-and two-dimensional results are essentially the same and 

R^sl.O (21) 

provided that b/2L « .56 an( j k » <1.0 


The wave number k^ 1 must be expressed in terms of airplane temporal 
frequency and is commonly related to the reduced frequency 


to.c 


k^’ = -^-.(particularly for calculations based on a one-dimensional 


turbulence model) as follows 


2L 


V = k i - = 

c 


( 22 ) 


As indicated in previous equations; e.g., equation (l6), the parameter 


■TT- is introduced by the two-dimensional turbulence model. Thus both chord c 

CU 


and span b are involved and, hence, aspect ratio should be explicitly 
identified, lest its value be inadvertently changed as variations in 


2^- are considered. Therefore, equation (22) is written as 


k * = k n ^ AR 
1 Id 


(23) 


For "k ' < 1.0 and ^ « .56, Equation (23) becomes 


(2k) 


k, AR < ~ « .56 



This indicates the range of k^ for which RCfcl.O. 

The relationship in equation (24) Vas refined by a numerical evaluation 
of the integral of I g. It was found that 


.95 < R < 1.0 

for values of b/2L up to 7t/20 or 0.l6. Equation (24), there.fqre, can 
be expressed as 

k 1 AR <,l6 

An example of the relationship between the exact integrand Ig, and 
the approximate integrand 1^ is given in figure 2, with linear scales, 
for k^ ' = .5 and b/2L = ir/20. It can be seen that the effect of S 
is very small for this case. 

High frequencies .- Consider now the frequency range k^ f » 1,0, 

From equation (20), the two-dimensional gust spectrum in this frequency 
range is 



A typical spectrum and the S function are shown individually in figure 3 
as a log-log plot. The shape of the spectrum as a function of kg* is 

dependent only on the value of k^ f . The break frequency is 


For 


Therefore 


^/ B r ' - 4 ,' AR ' 


y « y 

B,I 

T ^ *5 

h * ir Y 3 
3 1 0 =■ 


^kfr 3s 


The integral of S with respect to kg* converges rapidly, reaching 
90 percent of the maximum value at an upper limit of 
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It follows, then, that if 


X - 4 tt ar >* R7r ¥ 


A iR » 2 tt 


then T (X A ) a v/' 

-“V '^V TrJ 0 „4; 3 

The exact expression (from equations (13 ) and. ( 20 )) is 




JL fXf+x'l 5 M' 


The integrands cf equations (27) and (28) are plotted with linear scales in 

V 

figure 4 for k ’ = 100 and = — . It is seen tnax the differences are 

very small and so equation (27) is valid for large values of k^ 1 . 

The combination of equations (27) and (l6) (with a change in variable of 
integration) yields 


** * & ? r sj **- 




where 


r*«*- -rx j * ' * 


Prom equation (34), for k^ 1 = k^. AR » —, 

V „ 


iL ex) 


The weighting function R is therefore 

R ~2L jr. 

R ~ b 2k ' 
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R(k 1 ), 


TT 

2 AR k. 


to to 

provided that — < 2ir (This restriction on — is due to the initial 


k^ » 1,0 


k i AR » £ 


together with the requirement in equation (2fi). It is of no consequence 
to the present study, inasmuch as only values of b/2L less than sixteen 
percent are of interest ) . 

It is noteworthy that R in equation (2 9) does not depend on b/2L 
at all; it applies for all values of b/2L < 2ir. However, aspect ratio 
is an important parameter and R varies inversely with k^. 

Two-dimensional Gust Weighting Function .- The expression for R in 
equation (29) for the high frequency asymptotic condition (lc » 2 $) 

AR 

together with the results from consideration^^ th|e low frequency asymptotic 
condition (equation (2l)), R-^l.Q for k^ < — ^ , allows the construction 

of a function that approximates R over the entire range of k^. The 

following equation satisfies these asymptotes. 


*(■*;) 


I •/- £ A R. ! 
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This function is plotted on figure 5 and is identified by the label "Dryden, 
It is of interest to observe that equation (30) is completely independent of 
b/2L for values of b/2L less than 0,16, Equation (30) is applicable, 
therefore, to very small airplanes. The amplitude of R is nearly unity 
at low frequencies and is attenuated inversely with frequency at the high 
frequencies. The attenuation is governed, otherwise, only by aspect ratio 
as indicated by the break frequency of the asymptotic lines at 
. _ tt . The attenuation at a given frequency increases with an increase 


in aspect ratio. 

Equation (30) provides results similar to those.. obtained by numerical 
integration in reference 11 and by asymptotic analysis of cross spectra in 


reference 13. However, in both references, the results are described by 
families of curves rather than by a single function. In reference 11, 
aspect ratio is not explicitly identified, and in reference 13, a separate 
function for each aspect ratio is calculated. Thus, the fundamental role 
of aspect ratio is not described. Furthermore, although calculations in 
both references are made for values of b/2L as low as 0 . 05 , evidence is 
not presented to indicate that the attenuation at high frequencies exists 
for all small values of b/2L, including indefinitely small values. 

The similarity of results obtained by the method described herein to 
those in references 11 and 13 was verified by generating comparable families 
of curves from a weighting function based on the von Karman turbulence 
spectrum, used in the references. 


i + -=■ tf 

1.339* 1 (31) 

This function is also plotted on figure 5. It differs from equation (30) 
by a slightly lower break frequency. In spite of the fact that the weighting 
function is exact only for the very low and the very high frequencies, the 
results differ from those in reference 11 by a little more than 10 percent at 
the most ; the differences are even smaller with respect to results in 
reference 13. 

The weighting function R is useful primarily as a tool for revealing 
both the characteristic attenuation of lift at high frequencies for arbi- 
trarily small values of b/2L and the role of aspect ratio. VThile R may 
be suitable to obtain quick approximate results for conceptual design pur- 
poses (for unswept wings), more accurate results can be obtained by use of 
computer solutions for the aerodynamic forces on lifting surfaces of 
various shapes (e.g., kernel function, ref. 16 ; and doublet-lattice, ref. 17). 
The turbulence spectrum in terms of dual wave numbers, used herein, is 
particularly well suited for use in such calculations. 

The procedure used herein can be applied to tapered unswept wings, 
although more effort is required to obtain closed- form integrals. From 
qualitative considerations, the effect of taper is similar to a reduction 
in aspect ratio of an untapered wing (due to loss of effectiveness of the 
outboard portion of the surface ) . The result is an increase in the break 
frequency of the weighting function R, Quantitative evidence of this is 
found in reference 11. Triangular span loading (a taper ratio of zero) 
increases the break frequency by about 25 percent over that for constant 
span loading (taper ratio of one). This is equivalent to a 20-percent 
reduction in aspect ratio. 


NORMAL ACCELERATION SPECTRA 


The full significance of the ratio R is brought out by consideration 
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of the spectrum of the acceleration response of an airplane. This spectrum 
is described by combining equations (12), (13) » (l4), and ( 30 ) which yields 



• The weighting function R^) is shown in the above equations as a separate 
entity. It could be logically combined with the turbulence spectrum to 
define an effective spectrum for the two-dimensional turbulence model as in 
equation (l4). However, R is a function of airplane parameters only, notably 
aspect ratio and so, also logically, it could be combined with the airplane 
frequency response function. 


Functions The important properties of the airplane acceleration 
spectrum can be illustrated through the use of rather simple approximations 
for the several functions of frequency in equation (33). 

The Dryden spectrum for the one-dimension turbulence model is 



This is further simplified by use of asymptotic functions to 



The airplane frequency response function considered is for a rigid or 
quasi-flexible airplane with a single degree of freedom in plunging motion 
(Pitching motion is suppressed). The unsteady aerodynamic forces due to 
motion are ignored. The equation of motion is 
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The modulus-squared of the frequency response function can be written as 
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The frequency response function for this condition contains the significant 
features of the more general condition. 


The unsteady gust function selected is 
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This is an approximation of the unsteady gust force for a lifting surface 
having an infinite aspect ratio in incompressible flow, as used in reference 
3. Different approximations, but also based on infinite aspect ratio, were 
used in the studies of responses to two-dimensional turbulence in references 
10 and 11. This practice carries the tacit assumption that the variation 


of total lift with frequency 


is independent of the shape of the spanwise 


distribution of lift, which, of course, is a function of frequency 


This assumption is consistent with others inherent in a simple strip theory 
approximation for the aerodynamic forces, but apparently the validity has 
not been assessed. It would seem desirable to investigate the character- 
istics of unsteady gust forces by calculating them by means of the doublet- 
lattice or similar computer program. 


The four functions of frequency in equation (33) as described by 
equations (30) and (35) - (37) are illustrated in figure 6 as log-log plots 
of the asymptotic functions. The airplane acceleration spectrum is obtained 
from the sum of these logarithmic plots and is illustrated in figure 7 . The 
slopes of the logarithmic lines, the maximum amplitude, and the break-points 
of the lines are identified. 


The high frequency asymptote is of particular significance, and it is 
notable that for an aspect ratio of TT^('ilO) , the break-points of the R 
weighting function and the unsteady gust lift function coincide; thus, each 
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function contributes equally to the high frequency characteristic. The 
significance of the high frequency asymptote is associated with its influence 
on the second moment of the response spectrum that in turn influences the 
value of the average rate of zero crossings N as indicated by equation (2). 
The second moment of the acceleration spectrum°is shown as a log— log plot of 
asymptotic lines in figure 8. The corresponding function for the one- 
dimensional turbulence model (with B = 1,0) is shown by the dashed line. 

The integral for the one-dimensional turbulence model is infinite due 

to the asymptotic slope of minus one. Calculations of N for this case 

0 

have required truncation of the integration at an upper frequency limit 
corresponding to some stipulated percentage (close to 100$) of the mean- 
square value of the response as suggested in reference 2, In contrast, the 
integral of the function described by the solid lines is finite, thus 
offering the promise of determining N q without truncating the integration. 

This promise is evaluated by an application to a representative example. 


Application to Example Airplane 

The example airplane is closely akin to a short-haul transport in 
current service. The airplane features an unswept, nearly constant-chord 
wing. The pertinent aircraft and atmospheric parameters are listed in 
table 1 for the cruise condition. The pitching motion is suppressed. 

Note that the aspect ratio is 10 which results in about equal attenuation 
by the two-dimensional turbulence and the gust unsteady lift. 

A An “* 1116 s P ectra for normal acceleration at the center of gravity 

are presented in log-log plots together with the asymptotic lines in 
figure 9 and in linear plots in figure 10. The figures include results 
from the simplified method for the two-dimensional turbulence model and 
from the conventional one-dimensional model. The difference in associated 
values of A is about 8 percent. 

N. The spectra for the second moment are presented in like forms 

An 

in figures 11 and 12. The value of for the two-dimensional turbulence 

model is .68 per second which is believed to be a reasonable value. Of 
course, for the one-dimensional turbulence model the value of is 

infinite unless the integration of the second moment of the spectrum is 
truncated. 

Vanishing b/2L 


Turning new from the characteristics related to an airplane at a low but 
fixed value of b/2L» attention is given to the effects of b/2L approaching 
arbitrarily small values. Much useful information is obtained by recon- 
sidering the properties of the several frequency functions illustrated in 




figures 6-8. It is seen that the parameter b/2L appears only in the 
turbulence spectrum function. Consequently, it is convenient to combine 
the other functions into a single effective acceleration response function 
and to examine its relation to the turbulence spectrum for various values 
of b/2L. This iB illustrated in figure 13. The acceleration spectra at 
the bottom of the figure indicate that the maximum amplitude decreases 
directly with a decrease in b/2L and so approaches zero as b/2L approaches 
zero. Although the bandpass characteristics broaden toward the lower frequen- 
cies with decreasing values of b/2L, recall that this plot is logarithmic, 
and the asymptotic lower frequency is zero. The area under the acceleration 
spectrum, hence, the value of S’, must approach zero as b/2L approaches 
zero , provided that K and AR are held constant . 

In contrast to A, N 0 approaches a constant value as b/2L approaches 
zero. «_ is proportional to the radius of gyration of the response spectrum 
with respect to the zero frequency axis based on linear scales (See equation 
(2)). From the log-log plot on figure 13 it can be seen that the shape of 
the acceleration spectrum remains constant, except at low frequencies, as 
b/2L decreases. It follows that the radius of gyration (for linear scales), 
hence Nq will decrease somewhat and approach a constant value as b/2L 
becomes indefinitely small. 

It would appear trivial to compare the characteristics for the two- 
and one-dimensional turbulence models for b/2L - 0. However, there it' 
often some value in Interpreting asymptotic condition^ and this has been 
done in past studies; e.g., references 3 and 10, the latter being the moti- 
vation for the study presented here. As can be seen from the high-frequency 
asymptotes for the acceleration spectra in figure 13, the results for the 
one- and two-dimensional turbulence models are the same only for aspect ratio 
equal to zero, inasmuch as the break-point frequency associated with the R 
weighting function becomes infinite. The usual application of the one- 
dimensional model results in precisely this. The frequency lc is expressed 

• WlF x 

in terms of wing chord, i.e., k 1 = . At the same time the scale is 

stipulated to be some finite value while the span is assumed to be very 
small compared to scale and zero in the limit; the aspect ratio is thus 
forced to zero. 

Choice of parameters .- The use of the parameters b/2L and k^ in 

the final expressions (i.e., equations (30) and (35)-(37))is arbitrary. 

The parameter b/2L was selected because the study was directed toward 
effects for small values of the quantity. It appears in the results only 
in the expression for the one-dimensional turbulence spectrum as a function 
of k^. In view of the conclusion that the effects of the two-dimensional 

gust field are independent of b/2L for small values, the airplane size can 
be expressed in terms of a reference chord instead of span. The aspect 
ratio then would appear only in the weighting function R. 

It is felt that the choice of -k^ for the frequency expression produces 
simpler and more easily interpreted results than other alternatives such 
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as lii. The reduced frequency 1^ is a fundamental parameter widely used 

with regard to unsteady aerodynamic forpes, including those due to gust 
penetration. The frequency response function for the normal acceleration 
of the airplane is a simple expression in terms of k^. The shape of the 

function for a plunging airplane is dependent only on the mass parameter K. 
The weighting function R is dependent only on aspect ratio. As previously 
mentioned, the description of the one-dimensional turbulence spectrum in 
terms of is dependent only on airplane size and scale of turbulence . 

In general, if alternate expressions for frequency are utilized (e.g., ut), 
and if the airplane Bize is described by a reference chord, the weighting 
function R remains dependent on aspect ratio. 


CONCLUDING REMARKS 


The normal acceleration response of a rigid airplane with a constant- 
chord, unswept wing, flying in turbulence that varies in the spanwise and 
flight path directions, is analyzed. The analysis is restricted to a range 
of values of the ratio of wing span to turbulence scale from 30 percent to 
vanishingly small, a range which includes most airplanes. The analytical 
approach features an atmospheric turbulence spectrum in terms of wave numbers 
in both directions. The spanwise distribution of aerodynamic forces is 
approximated by strip theory. Spectral integrations with respect to spanwise 
wave numbers are obtained in closed form for asymptotic conditions attending 
small and large wave numbers in the direction of flight. 

Results show that the power spectrum of normal acceleration in the two- 
dimensional turbulence field is significantly attenuated at the higher 
frequencies relative to that for a one-dimensional field and that the attenua- 
tion persists for values of wing span to turbulence scale approaching zero. 

The attenuation is accounted for by a simple weighting function of frequency 
that is dependent only on the aspect ratio of the wing, the attenuation 
increasing With an increase in aspect ratio. 

The attenuation, together with that from the unsteady flow of gust 
penetration, permits the determination of the average number of zero crossings 
without recourse to an integral truncation that is often required in one- 
dimensional turbulence calculations. It is notable that for a wing aspect 
ratio of about 10, the contributions of the two dimensional gust field and 
gust penetration to the attenuation are equal. 

The above result is based on the use of an unsteady gust force for a 
lifting surface having an infinite aspect ratio, a frequently used 
approximation. This practice is based on the tacit assumption that the 
unsteady flow effects are independent of the shape of the spanwise lift 
distribution. The adequacy of this assumption needs to be assessed by 
application of unsteady lifting surface computational programs. 
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TABLE X.- AIRCRAFT AND ATMOSPHERE PARAMETERS 


Wing area, S 

39 m 2 . 

(420 square feet) 

Mean wing chord c 

1,98 m, 

(6,5 feet) 

Wing span 

19.8 m, 

(65 feet) 

Weight, W 

50,042 Newtons, (11,250 pounds) 

Altitude, h 

9XU m. 

( 3,000 feet) 

Cruise Speed, V 

80,5 m/s 

, (264 feet/second) 

Turbulence scale, L 

762 m, 

( 2,500 feet) 

Aspect ratio, AR, 

10 


Lift-curve slope, C^ , 

ot 

5 


Span-to-scale ratio, ~ 

.013 


Cruise lift coefficient., C* , : 

u n 

.353 


Mass parameter, K, 

94.1 
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